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Using a combination of first-principles total energies, a cluster expansion technique, and Monte 
Carlo simulations, we have studied the Li/Co ordering in LiCoC>2 and Li- vacancy/Co ordering in the 
□ C0O2. We find: (i) A ground state search of the space of substitutional cation configurations yields 
the CuPt structure as the lowest-energy state in the octahedral system LiCo02 (and DCoC^), in 
agreement with the experimentally observed phase, (ii) Finite temperature calculations predict that 
the solid-state order- disorder transitions for LiCoCh and □CoCh occur at temperatures (~5100 K 
and ~4400 K, respectively) much higher than melting, thus making these transitions experimentally 
inaccessible, (iii) The energy of the reaction Etot(o~, LiCo02) — E to t(o~, DCoC^) — £tot(Li, bcc) gives 
the average battery voltage V of a Li^CoC^/Li cell for the cathode in the structure a. Searching 
the space of configurations a for large average voltages, we find that <r=CuPt (a monolayer (111) 
superlattice) has a high voltage (V=3.78 V), but that this could be increased by cation randomization 
(^=3.99 V), partial disordering (V=3.86 V), or by forming a 2-layer Li2Co2C>4 superlattice along 
(111) (V=4.90 V). 



I. INTRODUCTION 

Much like the ABC2 semiconductors (A, B = Al,Ga,In 
and C = N,P,As, or Sb) exhibit nation ordering in 
a tetrahedrally-coordinated networktl, the LiAfC>2 ox- 
ides (M=3<i transition metal)em form a similar series 
of structures based on the octahedrally-coordinated net- 
work with anions (O) on one fee sublattice and cations 
(Li and M) on the other (Fig. [l]). Cation arrange- 
ments in isovalent (III-III-V) or heterovalent (I-III-VI) 
semiconductor alloys have been observedEJ in the disor- 
dered, CuAu-type (CA), CuPt-type (CP), and chalcopy- 
rite (CH) structures (bottom row of Fig. [l]), while, cation 
arrangements in the oxides have been observedau in the 
disordered, CP, CH, D4, and Y2 structures (top row of 
Fig. [j]). Ab-initio total energy calculationso have shown 
that in the tetrahedrally-coordinated III-V semiconduc- 
tor alloys, the CuPt structure is the least stable [due to 
the fact that it represents a stacking along the elastically 
hard (111) direction], while the chalcopyrite structure is 
most stable (it possesses both the lowest electrostatic and 
strain energies) . Similar studies have been performed for 
the octahedrally-coordinated networks oi the spin alloy 
Mn^S-Mn^S and the lead chalcogenides.Q In this paper, 
we examine the energetics and thermodynamics of cation 
ordering tendencies in the octahedral LiCo02 oxide, and 
compare to the tetrahedral semiconductor case which is 
well studied. The LiCo02 compound i&.u&ed-^«a±|Md£j 
material in rechargeable Li batteries .m#ffl yOE3£3y 
When Li is de-intercalated from the compound, it creates 
a vacancy (denoted □) that can be positioned in different 
lattice locations. Hence, we will examine not only (a) the 
Li/Co cation ordering (different sites for Li and Co) prop- 



erties of LiCo02 (.TLi=l), but also (b) the vacancy/Co or- 
dering (different sites for □ and Co) in DCoC^ (xLi=0). 
A third type of ordering in these materials, vacancy/Li 
ordering in Li T ni_ x Co02 (0 < x^ < 1), is not treated 
here. 

Our calculation proceeds in three steps: (1) Total en- 
ergy calculations: We calculate the T=0 total energy 
of a set of (not necessarily stable) ordered structures 
via the full potential, all-electepH linearized augmented 
plane wave method (LAPW)Bli3 with all atomic posi- 
tions fully relaxed via quantum mechanical forces. We 
then PjOa JJias^-eiiergies onto a (2) cluster expansion 
f(7£;.BBBBBO This expansion is a generalized Ising- 
like expression for the energy of an arbitrary substitu- 
tional cation arrangement. Once the coefficients of the 
expansion are known, the Ising-like expression may be 
easily evaluated for any cation configuration. Thus, one 
can calculate (via first-principles) the total energy of a 
few cation arrangements, but then effectively search the 
space of 2 N configurations (where N is typically < 10 4 ). 
Specifically, the cluster expansion may be used to search 
the entire configurational space for stable ground state 
structures, where one can obtain low energy, but oth- 
erwise unsuspected states (i.e., states which are not in- 
cluded in the set of calculated total energies). Having 
obtained such a general and computationally simple pa- 
rameterization of the configuration energy, we subject it 
to (3) Monte Carlo simulated annealingEl to extend first- 
principles calculations (at zero temperature) to finite- 
temperatures, thus obtaining order-disorder transition 
temperatures and thermodynamic functions. 

We find for LiCo02 the following: 

(a) A ground state search of the space of substitutional 
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FIG. 1. Cation arrangements in tetrahedral (semiconductor) and octahedral (oxide) networks. The black and white atoms 
represent the cations, while the grey atoms are the anions. Shown are the names of the cation configurations, the structure 
itself (bold lines indicate superlattice planes), and the equivalent superlattice of cations. Note: Some of the structures are 
shown as only a portion of the complete unit cell. 



cation configurations yields the CuPt structure as the 
ground state in the octahedral LiCoC>2 system, in agree- 
ment with the well-established experimentally observed 
phaseJa We find that this result holds even if the CuPt 
structure is not included in the set of energies used to fit 
the CE parameterization. The CuPt cation structure is 
the least stable bulk structure in tetrahedral ABC^. 

(b) Finite temperature calculations predict that the 
solid-state order-disorder transition for LiCo02 occurs 
at temperatures (~5100 K) much higher than melting, 
thus making this transition experimentally inaccessible. 
In contrast, order-disorder transitionsjn isovalent tetra- 
hedral ABC 2 systems are < 1000 KB The addition of 
Li vacancies lowers this transition to ^4400 K; but, this 
transition temperature is still too high to be observed. 
Thus, the finite temperature calculations demonstrate 
that the observed disordered (rocksalt) phase of LiCo02 
is not thermodynamically stable, but is only stabilized 
kinetically. 

(c) The energy of intercalation reaction 
£ t ot(cr,LiCo0 2 ) - Stot(o-) n Co0 2 ) - £ to t(Li,bcc) gives 
the average battery voltage Vrpf a LirCoC^/Li cell for 
the cathode in the structure erjj thus providing a means 
for prediction of battery intercalation voltages from first- 
principles energetics. Searching the space of configura- 
tions a for large average voltages, we find that tr=CuPt [a 
monolayer (111) superlattice] has a high voltage (V=3.78 
V), but that this could be increased by cation random- 
ization (V=3.99 V), partial disordering (V=3.86 V), or 
by forming a 2-layer Li 2 Co204 superlattice along (111) 
(F=4.90V). 

(d) Ordered cation arrangements in LiCo02 are sta- 



ble, similar to the keterovalent tetrahedral I-III-VI2 (e.g., 
CuInSe2) systemsE-a, but is the opposite situation from 
the isovalent tetrahedral III-V systems such as GaInP2, 
in which bulk ordered compounds are unstable. 

(e) The relative order of structural energies in the oc- 
tahedral LiCo02 system is quite different from the tetra- 
hedral cases: £?(CuPt) < E(CH) < E (CA) in both 
LiCo0 2 and DCoC^, compared with E(CH) < E{Ck) 
< E (CuPt), universally found in the lattice- mismatched 
tetrahedral systems (Fig. ||). 

II. METHODS OF CALCULATION 

the cluster expansion (CE) technique JlBEES!! which 
consists of an Ising-like expression in which each cation 
is associated with the site of an ideal lattice (fee, in this 
case), and the pseudo-spin variable Si is given the value 
+1(— 1) ii &n A(B) atom is assigned to site i. Within this 
description, the energy of any configuration a of cations 
can be written as:t3 

E CE (a)=Y,D f J f Ti f (a), (1) 
/ 

where / is a figure comprised of several lattice sites (pairs, 
triplets, etc.), D f is the number of figures per lattice site, 
Jf is the Ising-like interaction for the figure /, and 11/ 
is a function defined as a product over the figure / of 
the variables Si, averaged over all symmetry equivalent 
figures of lattice sites. This expression incorporates the 
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TABLE I. Lattice averaged spin products IT/ (a) (for a few figures /) of the cation configurations a shown in Fig. |l| 
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effects of atomic relaxation (indeed, one does not require 
that the cations sit precisely at the ideal lattice positions, 
but merely that there is a one-to-one correspondence be- 
tween lattice sites and atomic positions) t5 We determine 
{J/} by fitting -Ece(c) of N a structures to LDA total 
energies -Elda given the matrices {H/(cr)} for these 
structures. Table | gives the values of the lattice-averaged 
spin products 11/ (a) and degeneracies Df for the struc- 
tures in Fig. |l|. The values of 11/ (a) sometimes take on 
interesting degeneracies: For example, CA and CH differ 
only by 11/ (a) for figures beyond the. nearest neighbor, 
and the CuPt and D4 structures haved'Ej equivalent pair 
correlation functions for all pair separations. Also, all 
odd-body correlation functions are zero for any structure 
which possesses A — > B (or Si — > —Si) symmetry, such as 
all those shown in Fig. [y This means that CuPt and D4 
possess not only equivalent pair correlations, but three- 
body correlations as well. Thus, in terms of Eq. ([j]), all 
terms in the expansion which correspond to one-, two-, 
and three-body figures do not distinguish between CuPt 
and D4, and thus the first cluster correlation which can 
break the degeneracy between these two structures is in 
the four-body terms. In spite of the similarities between 
the CuPt and D4 structures, the former is a superlat- 
tice along the [111] direction and hence is rhombohedral, 
while the latter is not a supcrlattice and is cubic (see 
Fig. [j]). Thus, the CuPt structure has one extra struc- 
tural degree of freedom (namely a c/a ratio) that the D4 
structure does not haveEO 

We use N a = 8 configurations in the fitting proce- 
dure. These are shown in Fig. [|. The choice of end- 
point configurations requires some discussion. The nom- 
inal end point configurations, LiO and CoO in the NaCl 
structure, do not obey the octet rule, as LiO has 7 va- 
lence electrons/formula unit, while CoO has (in addition 
to its filled t2 g shell) 9 valence electrons/formula. As 
a result, these nominal structures have a very high en- 
ergy. In the 1:1 structures (LiO)„(CoO)„, an electron 



will move from each CoO unit to fill the hole in the LiO 
unit, thus creating normal octet bonds. These "charge 
compensated" end-point compounds (LiO)* and (CoO)* 
will have a lower energy than the nominal LiO and CoO. 
Our calculations thus consider only charge compensated 
structures. Using the procedure of Wei et alB3 in treating 
heterovalent alloys, the conventional, high-energy "end- 
point" compounds LiO+CoO are not included in the CE 
because they are not charge compensated. Our CE could 
be used to predict the energies of (LiO)* + (CoO)*, 
and we will see that this energy is indeed lower than 
that of nominal LiO + CoO. We only include the eight 
(LiO) n (CoO) n compounds shown in Fig. [I] in our fit. 
These cation orderings correspond to both the observed 
structures in the LiM0 2 series (CuPt, D4, Y2, and CH) 
as well as other cation arrangements not observed in this 
series (CuAu,W2,V2, and Z2). For substitutional order- 
ing problems, it is possible to enumerate all configura- 
tions up to a given unit cell sizeBa The set of cation 
configurations considered here includes all of the possi- 
ble equiatomic structures with unit cell size up to eight 
atoms. 

The set of twelve figures / retained in the expansion 
is the "empty" figure, the first ten neighbor pairs, and 
the nearest neighbor tetrahedron. In fitting the LDA 
total energies to the cluster expansion, we include a La- 
grange multiplier with the constraint that the pair in- 
teractions should be as smooth as possible in recipro- 
cal space. This technique (more fully explained in Ref. 
^) allows one to retain more figures in the expansion 
than total energies, and also requires the pair interac- 
tions to be as convergent as possible in real space. Al- 
though more sophisticated versions of the cluster expan- 
sion approaches are available when one requires extreme 
accuracy and has access to a large database of structural 
energies, we use the simple real-space expansion of Eq. 
([j]) with the Lagrange multiplier. The use of this simple 
form is predicated on the assumption that if we specialize 
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TABLE II. Predicted structural information of observed LiCoC>2 phases. Where available, experimental data and other 
calculated results are shown. CuPt - ABC and CuPt - AAA refer to CuPt configurations of cations (alternately stacked Li/Co 
or n/Go close-packed layers) in ABC... or AAA... type stackings, respectively. For □CoCh, CuPt - ABC and CuPt - AAA are 
isostructural with CdCl2 and Cdl2, respectively. Bulk moduli for LiCo02 (CuPt) and CIC0O2 (CuPt) were calculated (present 
work) to be 2.4 and 2.8 Mbar, respectively. 



Compound 


Method 


a (A) 


c (A) 


Li-0 (A) 


Co-O (A) 


Vea (A 3 ) 

v cq j 


LiCo0 2 (CuPt - ABC) 


Expt. a 
FLAPW (present work) 
Pseudopot. b 


2.82 
2.81 
2.93 


14.04 
13.60 
13.2 


2.07 
2.08 
2.10 


1.94 
1.90 
1.96 


32.23 

31.2 

32.71 


□C0O2 (CuPt - ABC) 


FLAPW (present work) 
Pseudopot. b 


2.78 
2.88 


12.13 
12.26 




1.85 
1.90 


26.9 
29.36 


LIC0O2 (CuPt - AAA) 


FLAPW (present work) 


2.79 


4.74 


2.11 


1.90 


32.0 


□C0O2 (CuPt - AAA) 


Expt. c 
FLAPW (present work) 


2.822 
2.80 


4.29 
4.01 




1.91 
1.85 


29.6 
27.1 


LIC0O2 (D4) 


Expt. d 
FLAPW (present work) 


8.002 
7.90 




2.06 
2.05 


1.95 
1.91 


32.0 
30.7 


□C0O2 (D4) 


FLAPW (present work) 


3.85 






1.85 


28.5 



a Ref. 13 
b Ref. ~ 
c Ref. il 
d Ref. 4C 



to fixed composition (e.g., either x-i,s=\ or £Li=0) the ex- 
pansion converges quickly with a small number of terms. 
For our CE, the error in the input energies used in the 
fit is a negligible amount, < 1 meV/formula unit. To 
obtain some idea of the errors involved in the CE pre- 
dictions, we have removed some structures and figures 
from the fitting process, and examined the resulting er- 
rors: Removing the CuPt (CH) structure and the four- 
body tetrahedron figure from the fit produces an 11 (49) 
meV/formula unit error in the energy of CuPt (CH), neg- 
ligible changes in the other fitted energies, and an energy 
of the random cation arrangement which changes by only 
1 (7) meV/formula unit. The magnitude of these errors 
is quite small in terms of the energetic scales of cation 
ordering (~1000-2000 meV/formula unit) and Li interca- 
lation (~4000 meV/formula unit). 

The expression of Eq. (|l|) can be applied to differ- 
ent ordering problems, with a separate expansion con- 
structed for each situation. Here, we construct three sep- 
arate cluster expansions to describe three different types 
of structural energetics: 

(a) Formation enthalpies for different Li/Co arrange- 
ments a on the fee lattice: 

AH f (a, LiCoO?) = E tot {a, LiCo0 2 ) - £tot(LiO, Bl) 
-£ to t(CoO,Bl) (2) 

where the last two terms refer to LiO and CoO in the 
NaCl (Bl) structure with lattice constants obtained by 
minimizing the respective total energies with respect to 



hydrostatic deformation. The resulting CE will reveal 
Li/Co ordering tendencies at Xu=l. 

(b) Formation enthalpies for different n/Co arrange- 
ments a on the fee lattice: 

AHfja, aCgp 2 ) = E iot (a, □ Co0 2 ) + E tot (Li, bec) 
-E tot (UO,Bl)-E tot (CoO,Bl) (3) 

where £ , to t(Li, bec) is the total energy of Li in the bec 
structure with lattice constant obtained from total en- 
ergy minimization. The resulting CE will reveal D/Co 
ordering tendencies at xu=0- 

(c) The Li battery intercalation reaction energy for dif- 
ferent Li/Co (and D/Co) arrangements on the fee lattice: 

Aff rea ct(^) - E tot (a, LiCo0 2 ) - E tot (*-, DCoOa) 

-£ tot (Li,bcc) (4) 

<AH rcact is the energy gained upon complete de- 
intercalation of Li from LiCo02, relative to Li metal, and 
is simply the difference between Eq. ^ and Eq. (||). If 
one assumes that the Li is removed without a change of 
the cathode structure a (a topotactic reaction), 

LiCo0 2 (cr) -> Lii-^CoO^cr) + a;Li+ + xe~, (5) 

then (see, e.g., Ref. ||) the reaction energy A7J rcact of 
Eq. (Q) is equal to the integral of the (zero temperature 
and pressure) open circuit voltage V of a Li,;Co0 2 /Li cell 
between Li compositions XLi=0 and xu=\; 
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FIG. 2. Formation energies of cation ordering in tetrahe- 
dral and octahedral ABC2 compounds in structures shown in 
Fig. [j]. For the GaInP2 and CdZnTe2 compounds, the en- 
ergy scale was multiplied by 5 for visual clarity. The GaInP2, 
CdZnTe2, and CuInSe2 energetics were taken from Refs. 
|3!| |24|, respectively. The LiCo02 and DCo02 energetics are 
from the present work. "PS" represents the energy of a 
phase-separated mixture of AC+BC rocksalt (zincblende) bi- 
naries in the octahedral (tetrahedral) systems. "Rand" is the 
energy of a phase in which cations are distributed randomly 
(i.e., with no correlations) on their sublattice. 



-F 



dxV(a,x) = -V(er) 



(6) 



where F is the Faraday constant. Hence, |AiJ rea et| is 
simply the intercalation voltage averaged over Li compo- 
sition. We note that the intercalation voltage calculated 
in this manner is a bulk, thermodynamic quantity and 
does not contain contributions from the cathode surface 
or from kinetic phenomena. 

The total energies needed for Eqs. (§)-(§) have bcea, 
obtained using the first-principles full-potential LAPWtS 
method. In the LAPW calculations, we used the ex- 
change correlation of Ceperlejf and Alder as parameter- 
ized by Perdew and Zunger.c3 LAPW sphere radii were 
chosen to be 2.0, 2.0, and 1.3 a.u. for Li, Co, and O, 
respectively. A well converged basis set was used, corre- 
sponding to an energy cutoff of 25.5 Ry (RK ma ^— 6.57). 
Tests were performed placing the Co ip levels in a sepa- 



rate semi-core energy window as opposed to treating the 
Co 3p as a core state; negligible differences were found, 
and thus the latter was used in all the calculations de- 
scribed below. Brillouin-zone integrations are performed 
using the equivalent k-point sampling method, using k- 
points for each structure corresponding to the same 28 
(6x6x6) special k-points for the fee structure.Eil All total 
energies are optimized with respect to volume as well as 
all cell-internal and -external coordinates. Convergence 
tests of the energy differences (with respect to basis func- 
tion cutoff, k-point sampling, and muffin-tin radii) indi- 
cate that the total energy differences are converged to 
within -0.01-0.02 eV/formula unit. 

Spin polarized calculations were performed for LiCo02 
and DCo02 in the CuPt cation arrangement in both fer- 
romagnetic (FM) and anti-ferromagnetic (AFM) geome- 
tries. For LiCo0 2 , both the FM and AFM calculations 
converged to the non-magnetic solution (/ico=0). How- 
ever, for CIC0O2, both FM and AFM calculations showed 
a weakly magnetic solution (nco ~0.45 for both FM and 
AFM) with the total energy of the FM (AFM) state being 
14 (—0) meV/formula unit below the non-magnetic state. 
Because spin polarization only has a small effect on the 
energy of these compounds, the calculations below for the 
energetics of cation ordering in LiCo02 and 00002 are 
non-magnetic (NM). NM, FM, and AFM (with the ob- 
served alternating [111] layers of spins) calculations were 
also performed for CoO, with the AEM solution being 
lowest in energy, and hence used hereEd 

Having obtained the coefficients {<//} of the CE of 
Eq. (|l]) for the three types (a)-(c) of ordering reac- 
tions, we subjected Ece{&) to a Monte Carlo simu- 
lated annealing method for treating the configurational 
thermodynamics.^ A system size of 16 3 = 4096 atoms 
(with periodic boundary conditions) was used in all cal- 
culations. Monte Carlo simulations were performed in 
the canonical ensemble with the transition temperatures 
being calculated from the discontinuity in the internal en- 
ergy as a function of temperature, and the ground states 
determined from the simulation at a temperature where 
all configurational changes proved to be energetically un- 
favorable. This gives: (i) the T=0 K ground state struc- 
tures (from a simulation of a finite size cell initially at 
high temperature, and subsequently slowly cooled to a 
low temperature where all configurational changes proved 
to be energetically unfavorable), (ii) the pair correlation 
functions or atomic short-range order present in the dis- 
ordered alloy, and (iii) the order-disorder transition tem- 
perature, T c . 



III. T = FORMATION ENERGIES 

A. Energetics of Li/Co ordering in LiCo02 

The formation energies [Eq. (|J)] of LiCo02 in various 
cation arrangements are given in Table HI and calculated 
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TABLE III. FLAPW calculated formation energies (eV/formula unit) of various cation arrangements in LiCo02 and 
□ Co02+Li(bcc): AHf(a, LiCo02), AHf(a, DCoC^), (formation energies of a) and A// re act(cr) (average intercalation volt- 
age of LiCo02 relative to Li) are denned in Eqs. (^-(Q). V cq is the equilibrium volume (A 3 /formula unit) of LiCo02, and 8V 
is the change in volume upon Li extraction [i.e., T4 q (LiCo02) - \4 q (DCo02)]. All energies of various ordered, disordered, and 
partially-ordered cation arrangements in LiCoC>2 are from cluster expansions (CE) of FLAPW energetics, and are described in 
the text. 
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2.2 
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CuPt(r ? =0.88,SRO=0) 


-3.22 


+0.60 


-3.82 






D4(t ? =0.88,SRO=0) 


-3.21 


+0.71 
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structural properties are shown in Table [y. We note that 
the D4 structure is only slightly higher in energy than the 
CuPt structure. This competition is interesting because 
LiCoC>2 has been synthesized in tbc-. fl4rf t r lf c f H e by 
solution growth at low temperature BQuaE^BaEa (Al- 
though there was initially some discussion in the litera- 
ture about this low-temperature synthesized phase being 
CuPt with imperfect long-range orderO, i^is Jieiu-estab- 
lished that this phase is D4 (or "D4-like" ) The 
near degeneracy of the calculated energies of the CuPt 
and D4 structures is simply a consequence of their identi- 
cal pair and three body correlations II f (a) noted above. 
We will see that the four-body interaction J4 which dis- 
tinguishes these structures is quite small, consistent with 
the small energy difference between CuPt and D4. 



B. Energetics of D/Co ordering in DCo02 

The formation energies [Eq. ®] of DCoCb m various 
□ /Co arrangements are also given in Table IIP These 
configurations correspond to various arrangements of Co 
and □. We note the following: 

(1) The relative order of energetics is similar in DCoC^ 
as in LiCo02. There is only one qualitative difference: 
CH drops in energy significantly upon extraction of Li, 
and is lower in energy than the Y2 structure, whereas the 
reverse is true for LiCo02- 

(2) The separation in energy between CuPt and D4 
increases in DCo02 compared to LiCo02, due to the 
symmetry of the phases: Upon extraction of Li in the 
rhombohedral CuPt structure, the c/a ratio decreases 
significantly, providing a significant source of energy low- 
ering for DCo02 - CuPt. D4, on the other hand, is not a 
layered supcrlattice in any direction and has cubic sym- 
metry. Hence, the cell parameters of DCoC^ (D4) can- 



not distort in any preferred direction, and consequently, 
□ C0O2 (D4) does not relax as much as CuPt. 

(3) The CuPt structure of CIC0O2 (isostructural 
with CdC^) has an ABC... stacking of the cation 
planes. However, xecent electrochemical measurements 
of Amatucci et a/O have succeeded in completely de- 
intercalating Li from LiCo02, forming a DCoC^ struc- 
ture which is isostructural with Cdl2, with the stacking 
of planes in an AAA... arrangement (see Fig. |3|) which we 
call "CuPt (AAA)" . These two structure are not related 
to one another by substitutional degrees of freedom, and 
thus are not describable by a single cluster expansion. 
To examine these non-substitutional degree of freedom, 
we have performed total energy calculations of DCoC^ 
in both the CuPt and CuPt (AAA) structures (CdlA 
Consistent with the observations of Amatucci et aZ.til, 
we find that the DCoC^ in the AAA stacking is lower 
in energy than the CuPt structure by ~0.05 eV/formula 
unit. 

(4) We find that LiCo0 2 in the CuPt (AAA) structure 
(Fig. ||) is higher in energy than the CuPt structure (with 
ABC stacking) by ~0.15 eV/formula unit, in agreement 
with the fact that the observed CuPt ground state in 
LiCoC>2 has ABC stacking. 



C. Effect of cation arrangement on average voltages 

Table JnJ gives the calculated reaction energies given 
in Eq. (^for each of the cation arrangements a studied 
here. The average voltages for all cation arrangements 
considered are in the ~4 V range. In particular, the av- 
erage voltage for LiCo02 in the CuPt structure (3.78 V) 
is up, reasonable agreement with measured values (4.0 
V)B'E3|iJU an( i pseudopotential calculations (3.75 V) 
Some configurations like CH, show a marked relaxation 
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FIG. 3. Different stacking arrangements of close packed 
cation planes in LiCo02 and DCoC^. The white, black and 
grey atoms represent Li, Co, and O, respectively. Clock- 
wise from the top left, the structures are: (i) The stable 
LiCo02 (CuPt) phase (equivalent to the phase shown in Fig. 
[j] but from a point of view which emphasizes the layered na- 
ture of the compound) with close packed cation planes in 
the ABC... fee stacking, (ii) The □Co0 2 (CuPt) structure 
(isostructural with CdCl2) formed from extracting Li from 
the LiCo02 structure, with close packed cation planes in an 
ABC... stacking, (hi) The observed CIC0O2 phase "CuPt 
(AAA)" (isostructural with Cdi2) which corresponds to an 
AAA... stacking of close packed cation planes, (iv) A hypo- 
thetical LiCo02 structure "CuPt (AAA)", formed from inser- 
tion of Li into the □Co02(Cdl2) phase. 



of the nCo02 phase, and hence show a significantly lower 
voltage than the other configurations. -Jjiuis, as also has 
been pointed out by previous authorst-3c3, we find that 
first-principles calculations can provide predictions inter- 
calation energies and hence, battery voltages. 

An interesting aspect of the effect of cation ordering on 
average voltage is that LiCoC>2 in the V2 structure has a 
much higher average voltage than CuPt. This increase in 
voltage is of interest because V2 is a (LiO)2(CoO)2 (111) 
superlattice, whereas CuPt is a (LiO)i(CoO)i (111) su- 
perlattice. If one exchanges every other pair of cations in 
the CuPt layered sequence, the V2 layering is obtained. 
Thus, V2 is just CuPt with anti-sites on two out of every 
four layers. This suggests that anti-site defects Lic and 
CoLi in the LiCo02 CuPt structure, while energetically 
very costly, should increase the voltage of the compound. 
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FIG. 4. Pair interaction energies Jf and J(k) in both real 
(left) and reciprocal space (right). Interactions are shown for 
the cluster expansions of energies of Li/Co ordering in LiCo02 
[(a) and (b)], energies of D/Co ordering in CIC0O2 [(c) and 
(d)], and average intercalation voltage in LiCo02 [(e) and 
(f)]. In real space, positive values of J/ indicate a preferred 
tendency for unlike atoms ("ordering") and negative values 
indicate a tendency for like atoms ("clustering"). The labels 
a, P, 7, and 8 indicate the §[110], 5 [111], 5 [201], and [000] 
points in reciprocal space, which correspond to the ordering 
waves of the Y2, (CuPt and D4), CH, and phase separated 
structures, respectively. 



IV. CLUSTER EXPANSIONS OF Li/Co AND 
□/Co ORDERING AND AH Icact 

We can now use the set of first-principles calculated 
energetics described in Section III to determine a set of 
interaction coefficients of the cluster expansion [Eq. ([!])] 
for cation ordering. We have constructed three cluster 
expansions of three different types of ordering: (i) Li/Co 
ordering in LiCo02, (ii) D/Co ordering in CIC0O2, and 
(iii) cation ordering effects on the average intercalation 
voltage (see below). 

The pair interactions J/ found from our CE fits are 
shown in Fig. both as real-space pairs J(|Ri — Rj|) 
and as the lattice Fourier transform in reciprocal space 
J(k). We see that the pair interactions in real space 
are decaying with distance quite quickly, indicating con- 
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vergence of the expansion. In reciprocal space, the pair 
interactions in Fig. ^| show some interesting properties: 
Minima in J(k) indicate wavevectors where composition 
waves are likely to form low energy structures. The four 
body interactions J4 found from our CE fits are much 
smaller than the pair interactions (e.g., J2) with the ra- 
tio J 4 / J 2 =0.004, 0.02, and 0.04 for the types of ordering 
(i)-(iii) above. 

For Li/Co ordering in LiCoC>2, the three local minima 
of J(k) are located at three wavevectors (shown by bold 
arrows in Fig. |): L-point 1(111), W-point ±(201), and 
near the K-point ^(110). (Additionally, an extremely 
shallow minima occurs between the T and X points.) 
These three wavevectors are the composition waves used 
to build all of the structures in the LiMC>2 series: CuPt 
and D4 [±(111)], CH [±(201)], and Y2 [±(110)]. The 
global minimum of J(k) occurs at the L-point [^(111)], 
the composition wave used to construct the CuPt and 
D4 structures. Thus, we anticipate that the pair interac- 
tions J(k) for cation ordering in other LIMO2 systems 
(for other transition metals M) is likely to be qualita- 
tively similar to that of Fig. |] with changes in the rel- 
ative minima at these three points. For D/Co ordering 
in nCo02, the minima in J(k) occur at the same points 
as in the case of Li/Co ordering in LiCo02, indicating 
that the relative ordering tendencies are similar in the 
two systems. 

For the cluster expansion of average intercalation en- 
ergy, the minimum of J(k) occurs at the T point, the ori- 
gin of reciprocal space (also shown by a bold arrow in Fig. 

5), indicating that phase separation into LiO+CoO 
should produce a low Ai? roact , and hence a high voltage. 

Once the interactions {J/} are obtained, Eq. (|]) pro- 
vides an efficient parameterization of the energy of any 
configuration. Applications of this cluster expansion pa- 
rameterization which we now discuss include a search of 
configuration space (via a Monte Carlo simulated anneal- 
ing algorithm) for ground state structures, which need 
not necessarily be included in the input set, thus open- 
ing the possibility of discovering unsuspected low energy 
states. One can also perform Monte Carlo simulations 
at finite temperatures to assess the thermodynamic and 
order-disorder properties of the system. Finally, one can 
easily calculate the energetics of disordered and partially 
ordered cation arrangements in these systems. 

As pointed out previously, the cluster expansions use 
as input only charge compensated compounds, and there- 
fore can be used to predict the energies of charge compen- 
sated (LiO)* + (CoO)*. We find from our CE of LiCo0 2 
that (LiO)* + (CoO)* is 0.79 eV/formula unit lower 
in energy than the nominal, non-charge-compensated 
LiO+CoO. Similarly, the CE of CIC0O2 predicts that 
(□O)* + (CoO)* is 0.84 eV below the non-compensated 
compounds. 



A. Ground States 

The simulated annealing algorithm finds the CuPt 
structure as the low temperature state. In Tabic III, 



we simply note that this structure was the lowest in en- 
ergy of the eight structures calculated by LAPW. But, 
the simulated annealing prediction of the ground state 
demonstrates that CuPt is also the lowest energy config- 
uration out of an astronomical number of possible con- 
figurations (without symmetry, there are ~ 2 N possible 
configurations that the algorithm could explore, where 
7V=4096). For our cluster expansion of nCo02, the sim- 
ulated annealing algorithm also finds CuPt as the lowest 
energy substitutional configuration. As we have already 
shown above, non-substitutional configurations are even 
lower in energy for the CIC0O2 system (e.g., the Cdl2 
structure). 

By combining the simulated annealing algorithm with 
the cluster expansion of average voltage, one can search 
for the cation configuration with maximum voltage. This 
search yields a phase separated (LiO+CoO) configura- 
tion (5.8 V). De-intercalating Li from this configuration 
would correspond to the artificial case of: LiO — > Li .i 
□ O(fcc). It is interesting to note that other authorsEEl 
by very different means have also arrived at the conclu- 
sion that this admittedly artificial case corresponds to 
the theoretical maximum voltage in LiM O2 compounds. 



B. Order-Disorder Transitions 

For LiCo02, the order-disorder transition between the 
low-temperature CuPt phase and the high-temperature 
disordered phase is predicted to occur at ~5100 K (Fig. 
||), well above the melting point of this material. (Note 
that the calculations in this paper are all solid-state, and 
thus do not consider the liquid phase.) Antaya et al.i3 
report a disordered rocksalt phase of LiCo02, grown by 
laser ablation deposition at 150° C, whereas growth at 
higher temperatures results in either the D4 or CuP-i 
phases. Our calculations indicate that the observedEil 
disordered rocksalt phase of LiCo02 is not thermodynam- 
ically stable, but is rather only stabilized kinetically, con- 
sistent with the fact that the disordered phase can only 
be grown at low temperatures. By performing a simula- 
tion of D/Co ordering DCo02 at finite temperatures, we 
were able to ascertain the effect of vacancies on the order- 
disorder transition temperature. Upon complete removal 
of Li, the order-disorder transition of nCo02 drops to 
^4400 K, still much too high to be experimentally acces- 
sible. Thus, the addition of Li vacancies is not likely to 
make the disordered rocksalt phase thermodynamically 
accessible to experiments (although this phase is still ki- 
netically accessible). 

At temperatures below 2000 K, the CuPt structure 
is predicted to be completely ordered. Thus, at any 
growth temperatures where thermodynamic equilibrium 
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FIG. 5. Monte Carlo calculated (a) energy relative to the 
T=0 energy of the CuPt structure and (b) heat capacity as 
functions of temperature for LiCo02- The transition between 
the CuPt and disordered (rocksalt) LiCo02 can be seen at 
~5100 K. 



is achieved, the CuPt phase should form with a long- 
range order (LRO) parameter of nearly unity. Thus, 
anti-site defects Lic or CoLi are probably not formed 
under conditions of thermodynamic equilibrium. Also, 
since CuPt is completely ordered by 2000K, even the 
D4 structure is not stabilized by thermodynamic fac- 
tors (i.e., thermal fluctuations in energy are smaller than 
the CuPt-D4 energy difference for temperatures of in- 
terest). However, the D4 structure has been observed in 
low-temperature solution grown and laser ablation-grown 
samples, which are probably not equilibrium phases. 



C. Properties of Disordered and Partially Ordered 
Cation Arrangements 

Using the CE, we can compute the energetics of any 
cation arrangement such as random alloys or any disor- 
dered (short-range or long-range ordered) phases. These 



energetics of several such phases in Tabic III 

Random alloy: The perfectly random alloy is a phase 
in which Li and Co atoms (or □ and Co for DCo02) are 
distributed on their cation fee sublattice with no atom- 
atom correlation between cation sites. This corresponds 
to the enthalpy as T — ► oo. The energy of this random 
phase is easily computed from the cluster expansion of 
Eq. ([!]), since the absence of atomic correlations leads 
to the simple values 11/ = 0, and thus the energy of 
the random alloy is given by -Ece (Random) = Jo. The 
energies of random cation arrangements in LiCo02 and 
□ C0O2 are shown in Table III. The ordering energy of an 
ordered compound a is the energy required to construct a 
from the random cation arrangement: 5E OI d(cr) — E(a) — 
-E(Random). From Table III, we can see that for both 
LiCo02 and DCo02, all ordered cation configurations 
considered have 5E ov d < 0, except for CA, V2, and Z2j_. 

Partial Short-Range Order. Because Antaya et a/.cll 
report the existence of some degree of CuPt-type or- 
dering in their disordered phase, we have also com- 
puted (Table III) the energetics of a disordered rock- 



salt phase with some degree of atomic short-range or- 
der (SRO). SRO is a finite-temperature effect, and is 
characterized in real space by the pair correlation func- 
tions IIo.„ + 1 for the nth atomic shell. Thus, the 
SRO parameters, IIo^i measure the extent to which spa- 
tial correlations exist in disordered alloys. The SRO 
parameters used to compute the energetics of for the 
first ten neighbor shells were obtained from a Monte 
Carlo simulation of the LiCo02 disordered alloy just 
above the order-disorder transition (in parentheses are 
the values for fully ordered CuPt or D4): -0.06(0.0), 
-0.27(-1.0), +0.03(0.0), +0.12(+1.0), +0.02(0.0), 
-0.07(-1.0), -0.02(0.0), +0.10(1.0), -0.01(0.0), and 
-0.01(0.0). Note that the energetic effect of SRO is 
to significantly lower the energy of the random phase in 
both LiCo0 2 and nCo02 by 0.27 and 0.40 eV/formula 
unit, respectively. 

Partial Long-Range Order: There have also been re- 
ports of long-range ordered LiCo02 (either CuPt or D4) 
with small quantities of Li on the Co sites, or vice versa. 
This amounts to a CuPt or D4 phase with partial long- 
range order (LRO). If the LRO parameter 77=1, then all 
Li and Co atoms reside completely on their own sublat- 
tice and LRO is perfect. However, for states of partial 
LRO, 77 < 1, and there is an amount of intermixing 
between sublattices. For simplicity, we assume that there 
are no short-range correlations between the intermixed 
atoms. In Table [II, we show the energetics of CuPt and 
D4 structures with LRO parameter r?=0.88, correspond- 
ing to 6% of Li on the Co sites, and vice versa. The 
LiCo02 energies of CuPt and D4 are both raised by 0.16 
eV/formula unit relative to the 77=1 fully ordered phases, 
while the corresponding increases for CIC0O2 is 0.20 and 







0.17 eV/formula unit. 

The cluster expansion of voltage can also be used to 
predict the average voltages of configurations not directly 
accessible to first-principles calculations (Table [II). In 
particular, we see that the random alloy (3.99 V) is pre- 
dicted to have a higher average voltage than the ordered 
CuPt phase (3.78-if). Since this phase has been produced 
by laser ablationpj it would be interesting to measure its 
electrochemical properties, in order to compare with our 
predictions. The increase in voltage due to disorder is 
significantly reduced when one considers the disordered 
phase with SRO described above (3.86 V). Thus, it is pos- 
sible that experimentally, the voltage of the disordered 
LiCo02 relative to CuPt could be used to indirectly de- 
termine the a mou nt of short-range order in the sample. 
Also in Table III are the voltages of partially long-range 
ordered CuPt and D4 phases with 6% of the Co atoms on 
the Li sites (77=0.88). Even this small amount of anti-site 
defects increases the voltages of CuPt and D4 by 0.05 and 
0.01 V, respectively. Note that for either LRO or SRO, 
the qualitative effect of disordering is the same: Disorder 
raises the energy of nCo02 more than LiCo02, and thus 
raises the average voltage. 



V. ENERGETICS OF OCTAHEDRAL VS. 
TETRAHEDRAL ABC 2 NETWORKS 

We now compare our results for cation ordering in the 
octahedral LiCo02 and DCoC^ systems with the well- 
studied cases of cation ordering in isovalent and het- 
erovalent tetrahedral ABC2 systems. In the (heterova- 
lent) octahedral LiM02 systems, ordered cation arrange- 
ments are stable with respect to LiO+MO. This ordered 
compound stability is qualitatively similar to heterova- 
lent tetrahedral systems, such as CuInSe2, where cation 
ordered phases are stable relative to decomposition into 
CuSe+InSe zincblende binaries. Just as in the LiCo02 
case, the binaries correspond to compounds (rocksalt LiO 
and CoO and zincblende CuSe and InSe) which do not 
satisfy the octet rule, and hence are relatively high en- 
ergy configurations. The high energy of the constituents 
"exposes" the ordered compounds (which do satisfy the 
octet rule) as stable in these I-III-VI2 systems. The sta- 
bility of ordered compounds in octahedral LiA/02 sys- 
tems is, however, in contrast to the isovalent tetrahe- 
dral A 111 B 111 CY systems, in which bulk ordered com- 
pounds are unstable with respect to AC+BC. In these 
III-III-V2 systems, the phase-separated state is the low 
energy ground state, and ordered compounds which have 
been observed have been shownEl to be a result of a com- 
bination of epitaxial strain and surface-reconstruction- 
induced ordering. 

The relative order of energies of ordered compounds in 
the octahedral LiCo02 and DCo02 systems is also quite 
different from the (isovalent or heterovalent) tetrahedral 
cases: E(CuPt) < E(CH) < E (CA) in both LiCo0 2 and 



□ C0O2, compared with E(CK) < E(CA) < E (CuPt), 
universally found in the lattice-mismatched tetrahedral 
systems. Also, in the octahedral systems, the ordering 
energy SE or d < for both CuPt and CH, while SE or d < 
for CH in the tetrahedral systems. 

The CuPt structure is preferred in octahedral networks 
due to strain energy arguments: For an octahedral ABC'2 
system which has distinct equilibrium A — C and B — C 
bond lengths, the CuPt structure has the property that 
the cell-internal distortion of the anions (C) accommo- 
dates any equilibrium A — C and B — C bond lengths, 
and maintains all A — C bond lengths equal to one an- 
other (and similarly for B — C). The D4 structure also 
possesses this optimal structural relaxation. This opti- 
mal bond-length accommodation is interesting in light 
of the fact that the cation CuPt structure in tetrahedral 
systems, when relaxed, possesses two equilibrium A — C 
bonds (and similarly for B — C) as opposed to the single 
A — C bond in the octahedral case. The distinction be- 
tween the two A — C bonds in tetrahedral CuPt is due to 
the fact that some of the A~C bonds in this structure are 
along the direction of cell-internal distortion and other 
A — C bonds are perpendicular to this direction. Thus, 
when C atoms are relaxed, these two types of A— C bonds 
adopt different bond lengths. However, in the octahedral 
CuPt structure, none of the A — C bonds are either along 
or perpendicular to the distortion direction of the anions, 
but rather all A — C bonds are at equivalent angles to this 
direction. Thus, when the anions relax, all A — C bonds 
are distorted by equal amounts. The calculated equilib- 
rium bond lengths in all cation arrangements for LiCo02 
are given m Table 0. One can see from this table that 
CuPt and D4 are the only structures for which there is 
only one type of symmetry-inequivalent Li-0 and Co- 
O bond. Structures with Li-0 and Co-0 bonds equal to 
one another (e.g., CH, CA) are energetically unfavorable. 
In tetrahedral systems, the configuration which possess 
the optimum structural geometry, analogous to CuPt in 
octahedral coordination, is the CH structure, which is 
the lowest energy ordered compound in size mismatched 
semiconductor alloys. This strain energy argument can 
also explain the relative stability of CuPt, CH, and CA 
in octahedral vs. tetrahedral systems: Using a simple 
valence force field which includes only energetic effects 
due to strain, one obtains the correct order of these three 
structures for both octahedral, and tetrahedral systems as 
compared with LAPW£3'E3a One should note, however, 
that in the LiM02 series, there are systems other than 
M=Co which possess ground states other than CuPt, 
e.g., the CH and Y2 structures. Thus, clearly strain-only 
arguments do not explain the totality of ordering tenden- 
cies in these compounds, as other effects must dominate 
in some systems. 

Another distinction between the ordering tendencies of 
the octahedral LiCo02 system with those of the tetrahe- 
dral systems is in the energy scale. In Fig. |[ the energy 
scale of the tetrahedral systems is multiplied by a factor 
of 5, and is still smaller than the octahedral energy scale. 
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TABLE IV. Calculated bonds lengths in various cation arrangements of LiCo02- 



Cation Structure Li-0 (A) Co-0 (A) 

CuPt 2.08 1.90 

D4 2.05 1.91 

Y2 1.93,2.08,2.28 1.86,1.92,1.94 

CH 1.94,2.24 1.89,1.94 

W2 1.93,1.97,2.09,2.12 1.86,1.90,1.93,1.98 

CuAu 1.87,1.99 1.87,1.99 

V2 1.89,2.21 1.84,2.04 

Z2 1.99,2.04,2.35 1.76,1.96,1.99 



The difference between the energy of the highest and low- 
est ordered compounds in the isovalent tetrahedral III- 
III- V 2 systems is &E(CuPt-CH)~0.1 cV/formula unit, in 
the heterovalent tetrahedral CuInSe 2 system it is SE(V2- 
CH)~0.7 eV/formula unit, whereas this difference in the 
octahedral systems is <lE(V2-CuPt)~1.4 eV/formula unit 
in the LiCo0 2 system and <5£(V2-CuPt)~2.4 eV/formula 
unit in CIC0O2. Thus the energetic effect of cation or- 
dering is much more dramatic in the octahedrally coor- 
dinated networks. 



VI. SUMMARY 

Using a combination of first-principles total energy cal- 
culations, a cluster expansion approach, and Monte Carlo 
simulated annealing, we have studied the cation order- 
ing in LiCo02 and □C0O2, and compared the ordering 
in these heterovalent octahedrally coordinated systems 
with previous studies of ordering in both isovalent and 
heterovalent tetrahedral ABC2 systems. We find many 
significant differences between ordering in octahedral and 
tetrahedral systems. In the heterovalent octahedral sys- 
tems, ordered compounds have negative formation ener- 
gies, and are hence stable. This is qualitatively similar 
to the heterovalent tetrahedral case, but distinct from 
the isovalent tetrahedral semiconductors, where ordered 
cation arrangements are unstable. Also, the relative or- 
der of ordered compound energies is different in the oc- 
tahedral systems studied here, relative to either isova- 
lent or heterovalent tetrahedral systems. In particular, 
for both the LiCo02 and DCoC^ systems, a simulated 
annealing ground state search of the entire cation con- 
figuration space yields the CuPt cation arrangement as 
the lowest energy ground state, whereas this structure is 
the highest energy configuration in tetrahedral III-III-V2 
systems. The scale of ordering energetics is dramatically 
different in the LiCo02 and CIC0O2 octahedral systems 
(<~ 1.5-2. 5 eV), compared with that of either heterovalent 
(e.g., ^0.7 eV in CuInSe 2 ) or isovalent tetrahedral semi- 
conductors (~0.1 eV). This difference in energy scales is 
also reflected in the different temperature scales of order- 
disorder problems in the two types of systems: While 
typical order-disorder temperatures in isovalent or het- 
erovalent tetrahedral systems are <1000 K, we find tran- 



sition temperatures of ^5100 K and ~4400 K for LiCo02 
and nCo02, respectively. 

Because LiCo0 2 is in a class of materials being stud- 
ied for use in rechargeable Li batteries, we have also 
examined the effects of cation ordering on Li intercala- 
tion energies and average voltages in Li^CoC^/Li cells. 
Searching the space of configurations a for large average 
voltages, we find that cr=CuPt [a monolayer (111) su- 
perlattice] has a high voltage (y=3.78 V), but that this 
could be increased by cation randomization (V—3.99 V), 
partial disordering (V=3.86 V), or by forming a 2-layer 
L12C02O4 superlatticc along (111) (17=4.90 V). 
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